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A toy model for the Blandford-Znajek mechanism is investigated: a Kerr 
^ , black hole with a toroidal electric current residing in a thin disk around the 

^v^ ' black hole. The toroidal electric current generates a poloidal magnetic field 
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field with remote charged particles, the rotation of the black hole and disk 
induces an electromotive force, which can power an astrophysical load at 
remote distance. The power of the black hole and disk is calculated. It is 
found that, for a wide range of parameters specifying the rotation of the 
black hole and the distribution of the electric current in the disk, the power 



/\ • of the disk exceeds the power of the black hole. The torque provided by the 



black hole and disk is also calculated. The torque of the disk is comparable 
to the torque of the black hole. As the disk loses its angular momentum, the 
mass of the disk gradually drifts towards the black hole and gets accreted. 
Ultimately the power comes from the gravitational binding energy between 
the disk and the black hole, as in the standard theory of accretion disk, instead 
of the rotational energy of the black hole. This suggests that the Blandford- 
Znajek mechanism may be less efficient in extracting energy from a rotating 
black hole with a thin disk. The limitations of our simple model and possible 



improvements deserved for future work are also discussed. 
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I. INTRODUCTION 

It's well believed that black holes exist in many astrophysical systems, such as in active 
galactic nuclei (AGNs), centers of galaxies, and some stellar binary systems. For a rotating 
black hole with an accretion disk, magnetic field threading the black hole and the disk could 
exist ( [|T| and references therein). Due to the rotation of the black hole relative to the 
magnetic field, the black hole's rotational energy can be extracted through the Blandford- 
Znajek mechanism [Q. Currents in the disk are needed to confine the magnetic field, but 
without using any mass accretion we can tap the huge rotational energy of the black hole. 
The Blandford-Znajek process provides a mechanism for extraction of energy from a rotating 
black hole which is more promising than the Penrose process in practice ^j^. For a long time 
the Blandford-Znajek mechanism has been considered as a reasonable process powering the 
radio jets in AGNs f^,^. Recently, this mechanism has been invoked in models for gamma- 
ray bursts (GRBs) where a rotating black hole with an accretion disk (or torus) forms 
through the collapse of a rotating massive star or the merger of a black hole with a neutron 
star [[?|-[T0[]. The Blandford-Znajek mechanism is favorable for these phenomena since ultra- 



relativistic Lorentz factors are required for jets in AGNs and GRBs and it's believed that 



very clean energy can be created through the Blandford-Znajek mechanism |11 



However, though some numerical calculations have been taken ( |T^ , |T^ and references 



therein), by now there are no clear answers to such questions as how the magnetic field is 
generated, how it is distributed, how much energy can be extracted from the black hole and 
disk, what portion of the extracted energy is attributable to the black hole, and if the energy 
extracted from the black hole is clean enough. Realistic cases are too complicated, we have 
very little observational clues. 

In this paper, instead of searching for complicated numerical solutions, we consider a 
semi-analytical toy model to probe the Blandford-Znajek mechanism. The toy model is a 
Kerr black hole with a toroidal electric current residing in a geometrically thin disk around 
the Kerr black hole. The toroidal current generates a poloidal magnetic field threading the 



black hole and disk. Due to the interaction between the magnetic field and remote charged 
particles, the rotation of the black hole and disk induces an electromotive force (EMF) on 
the black hole's horizon and on the disk. This EMF could be the energy source for remote 
astrophysical loads (such as the jets in AGNs and GRBs). The power and torque provided 
by the black hole and disk are calculated. It is found that for a wide range of parameters 
specifying the rotation of the black hole and the distribution of the current density in the 
disk, the power provided by the disk exceeds the power provided by the black hole. The 
torque provided by the disk is comparable to the torque provided by the black hole. This 
agrees with the results in |@,|14| and suggests that the Blandford-Znajek mechanism may be 
less efficient in extracting energy from a rotating black hole for the thin disk case. 

Though the case of thin disks is simple for calculation, in real astronomy accretion disks 
may be geometrically thick and the case of thick disks may be quite different form that 
of thin disks P, pT|p!5| . However the simple model presented in the paper gives a complete 
semi-analytical example: from the generation of magnetic fields to the extraction of energy 
from the black hole and the disk. The extension to the case of thick disks is extremely 
interesting and challenging. 

II. DESCRIPTION OF THE MODEL 

Realistic cases of a rotating black hole with a disk and magnetic field are likely to be 
extremely complicated. However, if the the magnetic field associated with the black hole 
has a somewhat poloidal structure, it can be modeled as being generated by some toroidal 
electric current outside the black hole's horizon. This toroidal electric current most likely 
resides in the disk around the black hole. Thus here we consider a model of a Kerr black 
hole with a thin disk in the equatorial plane and there is a distribution of electric current in 
the disk. The Kerr black hole has mass M, angular momentum Ma [throughout the paper 
we use the geometric units with G = c = 1 and the Boyer-Lindquist coordinates (t, r, 9, (p) 
for Kerr black hole]. Then the angular velocity of the black hole's horizon is 
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where rn = M + v'M^"^^^ is the radius of the Kerr black hole's outer horizon (a^ < M^). 
As usual the disk's angular velocity is taken to be the relativistic Keplerian angular velocity 
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The outer edge of the disk is at r = r^, the inner edge of the disk is taken to be at the 



innermost stable circular orbit (the "marginally stable" orbit) in the equatorial plane |[T6[| : 



m{3 + Z2- [(3 - zi)(3 + zi + 2z2)f^} , 



(3) 



where ^i = 1 + (1 - ayW)^/^ [(1 + a/My/"^ + (1 - a/My/"^] and Z2 = {'ia'/M^ + z^f'^. 
(The assumption of a Keplerian disk is valid only if the magnetic field is weak enough 
p!7| , [T8[| .) The toroidal electric current in the disk has a surface density J = J{r) which is 



distributed between r^s and r^ > r^s- 

The magnetic field generated by a single toroidal electric current at a fixed radius [i.e. 
the current density J{r) is a delta function] has been well investigated by many authors 
|T9|^T1. By linear superposition, the magnetic flux through a surface bounded by a circle 
with r = const and 6 = const is 



^(r, e) = 2nA4r, 6) = 27r J{r')-^dr', 
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where A^ is the toroidal component of the electric vector potential, and dA^/dr' is 
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where A = r^ - 2Mr + a^, S = r^ + 02005^6', A = (r^ + a^)^ - Aa^sin^6', m = (r 
M)/\/ M'^ — a?, Pi{z) and Qi{z) are Legendre functions, and Pi{z) = dPi{z)/dz, Q'i{z) 
dQi{z)/dz, and the coefficients a[, a^, /3[, and /?/ are respectively 
(1) for r < r', p[ = Pj = for all /; but 
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where A' = A(r = r'), S' = S(r = r', ^ = 7r/2) = r'^ and A' = A{r = r',9 = n/2). The 
normal component of the magnetic field in the disk is 

1 fA^'/' d^{r,7c/2)^ ^^^^ 
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The poloidal magnetic field on the horizon (which is perpendicular to the horizon) is 

1 d^{rH,e) 



B 



H 



(11) 



27r(r|^ + a2)sin^ dO 
Due to the interaction of the magnetic field with remote charged particles, the rotation 



of the black hole and disk induces EMFs on the horizon and in the disk [^ . The total EMF 
on the black hole is 
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£h = r^n^*(r^), (12) 
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where '^{th) = ^(^h,71"/2) is the magnetic flux through the northern hemi-sphere of black 
hole's horizon. The total EMF in the disk is 
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where \l/(r) = \E'(r, 7r/2), r^ < r^ is the radius within which the disk's energy is available. 
(As in the case of the Sun, the open magnetic field lines may emerge from only a small 



fraction of the disk surface p3|. So we should choose Vc -C r^ in practice.) These EMFs 
could be the energy source powering a remote astrophysical load. The black hole, the disk, 
and the remote astrophysical load form an electric circuit in series. (Alternative type of 
circuits are possible, see the discussions in Sec. |^) The circuit has two batteries — one is 
the Kerr black hole with EMF E^ and internal resistance Z^ (which is of several hundred 
ohms), the other is the disk with EMF Ed and negligible internal resistance (i.e. the disk's 
resistance is supposed to be <^ Z^) — and a resistor which is the astrophysical load with 
resistance Za = aZn- Suppose Sd and Sh generate a (single) poloidal electric current 
I = {Sh + Sd) I {Zh + Za), which flows along the circuit from the disk to the black hole, 
from the black hole up to the remote astrophysical load along the symmetry axis, from the 
remote astrophysical load down to the disk at the circle with radius Vc {vms < Tc < Tb) in the 
disk. The power provided to the remote astrophysical load, P = I^Za, is sensitive to the 
ratio a = Za/Zh- P reaches its maximum at a = 1. a = 1 is called the impedance matching 
condition [^. Define ^ = Sd/Sh, the ratio of the power of the disk Po = —dEo/dt = ISd 
to the effective power of the black hole Ph = —dEn/dt = ISh — I^Zh (where Ej:, is the 
energy of the disk. Eh is the energy of the black hole; Pd + Ph = P) is 

Pd ^ jl + a)^ .^^. 

Ph a-i ' ^ ' 

If ^ > a, the power of the black hole is negative or zero (while the power of the disk is always 
positive). If ^cr < ^ < «, where ^cr = '^/(2 + a), the power of the black hole is positive but 



less than the power of the disk, i.e. < Ph < Pd- If ^ < ^cr, the power of the black hole 
is bigger than the power of the disk. Thus, the parameter ^ gives a sensible measure of the 
relative importance of the disk's power and the black hole's power. If ^ > ^cr, the disk's 
power dominates. From the definition of .^cr, ■Ccr < 1 always. (For the case of P = Pm i-e. 
a = 1, we have ^cr = 1/3-) Therefore, if ^ > 1, we can conclude that the disk's power always 
dominates the black hole's power. 

The torque produced by the black hole is 

where Lh = Ma is the angular momentum of the black hole. The torque produced by the 
disk is 

Td = -^ = ^ [vP (rj - vp (r^,)] , (16) 

where L£, is the angular momentum of the disk, \l/(rc) — ^E'(rins) is the magnetic flux through 
the disk inside the circle with r = r^- The ratio of Td to Th is 

III. SOLUTIONS OF THE MODEL 

If we do not care about the absolute values of the power and torque of the black hole 
and disk, the problem is "self-similar" in the sense that the mass of the black hole and the 
absolute magnitude of the magnetic field are not important for us. Then, we are left with 
three adjustable parameters: (1) a/M, where M is the mass of the black hole and Ma is the 
angular momentum of the black hole; (2) the shape of the surface current function J = J{r)] 
and (3) r^/M, where r^ is the radius within which the disk's energy is available. We take 
J = Jo{r/rms)~"', where n is a positive dimensionless number, Jq is the magnitude of J at 
the inner edge of the disk. The radius Vc is usually thought to be several times larger than 



rff 0] but much smaller than r^ [^. Thus we take r^ = r^s + i^t"h, where k is a positive 



constant which has magnitude ~ 1. Then the three adjustable dimensionless parameters 
are a/M, n, and k (the value of Jq is unimportant for us). Clearly, if the magnetic field 
in the disk, B^, keeps the same sign everywhere (no reversal of magnetic field), the larger 
the radius Vc is, larger the power of the disk is, and thus less efficient the Blandford-Znajek 
mechanism is. To avoid overestimating the power of the disk by choosing a large r^, in our 
numerical solutions we take k, = 0.2. Even for such a small value of k, we will see that the 
disk's power still dominates the black hole's power. Thus, our numerical calculation is taken 
for K, = 0.2 and a series of values of a/M from 0.1 to 0.99, n from 1 to 3 (For a stationary 
Newtonian disk we have n = 1 P3). (For n > 1, the results are insensitive to the value of 
Th/M for sufficiently large Vf,.) 



During the numerical calculation, the summation Y,i over Legendre functions ||T9 
truncated at /max = 10. [For a realistic disk with thickness A^, the summation over / should 
be truncated at /max ~ (^^)^^-] The results are shown in Fig. |I] - Fig. |[ The results are 
insensitive to n. Fig. |I] shows the dimensionless parameter ^ = Sd/^h, which measures 
the relative importance of the power of the disk to the power of the black hole. For the 
impedance matching case with Za = Zh, the power of the disk always dominates the power 
of the black hole. For example, for a/M = 0.9 and n = 2, we have ^ = Sd/^h = 0.68, 
Pd/Ph = 4.3 [Eq. (|lD with a = 1]. Fig. | shows the ratio Pd/Ph for the case with 
n = 2 and Za = Zh- For a/M < 0.54, the black hole's effective power P^ is negative (i.e. 
PZh > Sh', Pd is always positive since we have assumed Zo = 0); for a/M > 0.54, both 
Ph and Pn are positive but P^ > Ph- Thus the disk's power dominates the black hole's 
power. Fig. ^ shows the dimensionless parameter rj = T^/Th, which measures the relative 
importance of the torque of the disk to the torque of the black hole. The torque of the disk 
is comparable to the torque of the black hole. Fig. ^ shows the normal components of the 
poloidal magnetic field on the horizon and disk, which use Bq = Jq/c as unit. The magnetic 
field is most strong at the inner edge of the disk. 



IV. DISCUSSIONS AND CONCLUSIONS 

From the simple model presented in the paper we have got some insight into the 
Blandford-Znajek mechanism. Since the electric current generating the magnetic field must 
reside outside the black hole and most likely resides in a disk around the black hole, as 
described in our thin disk model, the strength of the magnetic field at the inner edge of the 
disk should be bigger than that on the horizon |TB,n|. The solutions of our model show that 



this is true. Though the average strength of the magnetic field in the whole disk could be 
smaller than that on the horizon, the local strength of the magnetic field at the inner edge 
of the disk could be much stronger. In fact, in our solutions B^ir = rms)l Bh{0 = 0) has 
never got smaller than 1. Our solutions show that, for a wide range of parameters specifying 
the model considered in the paper, the power of the disk is stronger than the power of the 
black hole. The torque provided by the disk is comparable to the torque provided by the 
black hole. These results are physically plausible since (1) the electric current generating 
the magnetic field resides in the disk and thus the magnetic field at the inner edge of the 
disk (where the current density is most strong) should be stronger than that on the black 
hole's horizon; (2) rapidly rotating black holes exclude stationary and axisymmetric mag- 
netic fields |2^, as seen from Fig. ^ (left diagram); (3) the disk has a surface area larger 
than that of the horizon and thus the magnetic fiux through the disk could be larger than 
that through the black hole; (4) the disk is a perfect conductor with a negligible resistance 
while the black hole's horizon has a significant internal resistance of several hundred ohms 
which consumes the power of the black hole and disk. 

For the model considered in the paper the angular momentum and energy lost from the 
disk are important. As the disk loses its angular momentum, the mass of the disk gradually 
drifts towards the black hole and gets accreted. Ultimately the power at the remote load 
comes from the gravitational binding energy between the disk and the black hole, as in the 
standard theory of accretion disk, instead of the rotational energy of the black hole. This 
implies that the Blandford-Znajek mechanism is not likely to be efficient in extracting energy 



10 



from a rotating black hole with a geometrically thin disk and magnetic field. 

Though the simple model is simple and convenient for calculations and we believe that 
some essential features for the Blandford-Znajek mechanism have been included, the limi- 
tations must also be emphasized: (1) We have assumed that the disk is geometrically thin. 
But in real cases this may not be true, especially for subcritical accretion. There could 
be an alternative type of models where, in the limiting case, poloidal fiux is confined by 
and excluded from a funnel formed by an ion torus with ADAF or ADIOS solution where 
the disk's power is manifestly zero P,[Tl],^. (2) We have neglected the magnetic coupling 
between the black hole and the disk. In fact there could be strong coupling between the disk 
and the hole with the former acquiring energy from the work done by the later p], p6| , p7 



which might lead that some of the power of the disk effectively comes from the rotation of 



the black hole [0. In particular, loops of closed lines connecting the disk with the hole 
could exist and play important roles in transportation of angular momentum and energy. 
For a fast rotating black hole [ajM > 0.36) with a relativistic Keplerian disk, the rotation 
of the disk is slower than the rotation of the black hole. So the closed magnetic field lines 
will transfer angular momentum from the hole to the disk and reverse the accretion fiow 



. Accompanying the angular momentum, energy will also be transfered from the black 
hole to the disk via Poynting fiux. Then for the remote load the energy directly extracted 
from the black hole is decreased, but it may be compensated by the energy extracted from 
the disk. We can imagine that with suitable conditions a steady state could exist when the 
power from the disk balances the energy fiux from the black hole to the disk and then all the 
power of the disk effectively comes from the spin of the black hole. (3) We have supposed 
that the poloidal magnetic field is generated purely by the toroidal currents in the disk. By 
doing so we have neglected the effect of the currents induced by the charges in the magne- 
tosphere. This allows to calculate the poloidal magnetic field from the potential instead of 
solving the complicated Grad-Shafronov equation. With this approach the back-reaction of 
the induced charges and currents cannot be taken into account, but it is particularly useful 
for semi-analytical investigations. (4) We have modeled the toroidal magnetic field with 
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a global poloidal current loop flowing through the black hole, the accretion disk, and the 
remote load. This allows us to conveniently compare the power of the black hole and the 
power of the disk since they have the same unique current. The realistic case will be more 
complicated since poloidal currents flow into and out of the black hole and the disk diversely 
at all radius. An alternative type of current loops is that the poloidal current associated with 
black hole and the poloidal current associated with the disk flow separately, i.e. we have 
two distinct poloidal current loops. In this case the black hole and the disk have different 
loads and it's hard to compare their powers. However if the disk and the black hole have 
comparable EMFs and the impedance matching conditions are satisfled for both the hole's 
circuit and the disk's circuit, the power of the disk will exceed the power of the black hole 
since the disk has much smaller internal resistance. But in this case, even the black hole's 
power is only a small faction of the total power, it can still be important in practice if the 



energy from the black hole is very clean and not mixed with the energy from the disk |]TT 

In conclusions, for the simple model considered in the paper we have shown that the power 
of the disk dominates over the power of the black hole. This suggests that the Blandford- 
Znajek mechanism may be less efficient for extraction of energy from black hole with a thin 
disk. However, due to many simplifications we have made for the model, improvements 
are required for getting more thorough understanding of the Blandford-Znajek mechanism. 
Especially the generalization to the case of thick disk deserves considerations. 
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FIGURES 
FIG. 1. The variation of the ratio of disk's EMF, ££), to the black hole's EMF, E-h, with respect 

to the parameter a/M. The curves are plotted from a/M = 0.1 to a/M = 0.99 for different values 

of n: 1,2, and 3 [see the labels at the left end of each curve; the surface current density in the 

disk is J oc r^"]. The long dashed line represents EdJEu = 1, the short dashed line represents 

^d/^h = 1/3 which is the critical value for the most efficient case with Za = Zh (where Zh is the 

electric resistance of the black hole, Za is the electric resistance of the remote astrophysical load). 

The parameter EnjEu measures the relative importance of the power of the disk to the power of 

the black hole [see Eq. (pl|) in the text and the discussions below it]. If Ed/£h > 1, the power 

of the disk always dominates the power of the black hole. If £-d/£-h > 1/3, the power of the disk 

dominates the power of the black hole in the most efficient case with Za = Zh- 

FIG. 2. The ratio of the disk's power Pd to the black hole's power Ph for the case with 
n = 2 and Za = Zh- The vertical dashed line shows the position (where a/M = 0.54) for 
Ph = I£h - I^Zh = 0. For a/M < 0.54, the black hole's power Ph is negative; for a/M > 0.54, 
the black hole's power Ph is positive but less than the disk's power Pq. The shaded region 
(0 < Pd/Ph < 1) shows the case when the Blandford-Znajek mechanism is efficient, i.e. when the 
black hole's power dominates the disk's power. In the unshaded region, the disk's power dominates 
the black hole's power, the Blandford-Znajek mechanism is less efficient. None of the models is in 
the shaded zone, i.e. the B-Z mechanism is less efficient in all models considered in the paper. 

FIG. 3. The variation of the ratio of the disk's torque, Tq, to the black hole's torque, Th, 
with respect to the parameter a/M. The curves are plotted from a/M = 0.1 to a/M = 0.99 
for different values of n: 1,2, and 3 (see the labels at the end of each curve). The dashed line 
represents Td/Th = 1. If Td/Th > 1, the torque of the disk dominates the torque of the black 
hole; if T^/Th < 1, the torque of the black hole dominates the torque of the disk. 
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FIG. 4. The left diagram shows the distribution of the magnetic field Bh on the black hole's 
horizon for n = 1 and different values of a/M. The right diagram shows the distribution of the 
normal component of the magnetic field B£) on the disk for n = 1 and different values of a/M, 
drawn from r = rms{o) to r = 8M. The inner edge of the disk is at r = rms{0')^ as denoted with 
the thick dots at the left end of each curve. 
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